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ABSTRACT 

This paper is part of a series devoted to the investigation of a large sample of brightest 
cluster galaxies (BCGs), their properties and the relationships between these and the 
properties of the host clusters. In this paper, we compare the stellar population prop- 
erties derived from high signal-to-noise, optical long-slit spectra with the GALEX 
ultraviolet (UV) colour measurements for 36 nearby BCGs to understand the diver- 
sity in the most rapidly evolving feature in old stellar systems, the UV-upturn. We 
investigate: (1) the possible differences between the UV-upturn of BCGs and those 
of a control sample of ordinary ellipticals in the same mass range, as well as possible 
correlations between the UV-upturn and other general properties of the galaxies; (2) 
possible correlations between the UV-upturn and the properties of the host clusters; 
(3) recently proposed scenarios where helium-sedimentation in the cluster centre can 
produce an enhanced UV-upturn. We find systematic differences between the UV- 
colours of BCGs and ordinary ellipticals, but we do not find correlations between 
these colours and the properties of the host clusters. Furthermore, the observations do 
not support the predictions made by the helium-sedimentation model as an enhancer 
of the UV-upturn. 

Key words: galaxies: evolution - galaxies: elliptical and lenticular, cD - galaxies: 
stellar content - ultraviolet: galaxies 
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1 INTRODUCTION 

The UV-upturn, first discovered by Code & Welch (1979), 
is the rising flux with decreasing wavelength from 2500 A to 
the Lyman limit observed in the spectrum of some massive 
elliptical galaxies (see O'Connell 1999 for a review). 

In spiral and irregular galaxies, flux emitted in the UV 
acts as an excellent measure of the current star formation 
rate (Gil de Paz et al. 2007). In the case of (mostly qui- 
escent) early-type galaxies, the flux at wavelengths shorter 
than 1800 A (FUV) is believed to be produced by hot, low 
mass Horizontal Branch (HB) stars, although the mecha- 
nism with which an old population can develop hot stars 
is still hotly debated (see Yi 2008). For the non- quiescent, 
low redshift elliptical galaxies, analysis of UV wavelengths 
longer than 1800 A (NUV) is a useful indicator of residual 
star formation (Pipino et al. 2009). Possible solutions for 
the origin of the increase in FUV flux in quiescent early- 
type galaxies include metal rich stars that lose mass and 
become HB stars, metal poor stars that evolve to very high 
temperatures and onto the HB, or mass transfer in binary 
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systems that can expose the hot core of a HB star (Donahue 
et al. 2010). These theoretical scenarios make different pre- 
dictions about the correlations between the strength of the 
UV-upturn and the metallicity, as well as about the evolu- 
tion of this feature with redshift in passively evolving sys- 
tems. Observations of the evolution of the FUV-V colours 
of Brightest Cluster Galaxies (BCGs) with redshift seem to 
favour the single stars hypothesis above the binary theory 
(Ree et al. 2007). However, the metal poor and metal rich 
hypotheses also have problems. For example, a large fraction 
of metal poor stars are needed to explain the strength of the 
UV-upturn observed in elliptical galaxies and the metal rich 
scenario seems to require mass losses from Red Giant Branch 
stars much larger than those predicted theoretically. 

A popular third possibility, that seems to solve some 
of the problems associated with the HB hypotheses, is the 
enhanced helium abundance hypothesis. This hypothesis is 
motivated by globular cluster observations containing pe- 
culiar features in the colour magnitude diagram that can 
only be explained by the presence of super-helium-rich sub- 
populations (D'Antonafe Caloi 2004; Norris 2004; D'Antona 
et al. 2005; Lee et al. 2005; Piotto et al. 2005). These popu- 
lations could be a major source of FUV flux in quiescent el- 
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liptical galaxies, although studies suggest that the deduced 
value of the helium abundance is too small to be influen- 
tial on galactic scales (Lee et al. 2005; Yi 2008). However, 
in the centres of clusters, an extra source of helium could 
come from sedimentation processes (Peng & Nagai 2009). If 
the helium sedimentation mechanism is important, the UV- 
upturn strength in the BCGs could be greatly enhanced. 
In fact, it is known from observations of nearby early- type 
galaxy samples, that BCGs are most likely the strongest UV 
emitters among the quiescent early-type galaxies in each red- 
shift bin (Burstein et al. 1988; Boselli et al. 2005; Donas et 
al. 2007). 

If the mechanism producing the UV-upturn in BCGs 
is different than that of ordinary elliptical galaxies, then 
the results on the evolution of the UV-colours with redshift 
mentioned above will not necessarily be representative for 
the whole population of objects hosting old stellar popula- 
tions. In fact, Atlee, Assef & Kochanek (2009) obtained very 
different results for UV-colour evolution using a sample of 
ordinary bright elliptical galaxies. 

In this work, we want to explore the mechanism(s) pro- 
ducing the UV-upturn in BCGs and compare it with that 
producing the UV-upturn in elliptical galaxies by study- 
ing the relation between the UV-upturn strength and other 
properties of the galaxies. Several authors have found cor- 
relations between the strength of the UV-upturn and the 
Mg-index, an indicator of metallicity, in early-type galaxies 
(Burstein et al. 1988; Donas et al. 2007; but also see Rich 
et al. 2005). The question remains whether correlations be- 
tween the UV-upturn and the spectral properties of BCGs 
exist, which will provide us with further clues about the 
origin of the UV-upturn. 

This paper is the third in a series of papers investigat- 
ing an overall sample of 49 BCGs in the nearby Universe 
for which we have obtained high signal-to-noise ratio, long- 
slit spectra on the Gemini and WHT telescopes. This large, 
spatially-resolved, spectroscopic sample of BCGs allows pos- 
sible connections between the kinematical, dynamical and 
stellar population properties to be studied. The first paper 
was devoted to the spatially resolved kinematics of the BCGs 
(Loubser et al. 2008, hereafter Paper 1). In the second pa- 
per, we derived Single Stellar Population (SSP)-equivalent 
ages, metallicities and a-abundance ratios in the centres of 
the galaxies, and then systematically compared these with 
a sample of ordinary elliptical galaxies, as well as with the 
host cluster properties (Loubser et al. 2009, hereafter Paper 
2). 

Here, we turn our attention to the UV-upturn and its 
correlations with the properties derived from the optical 
spectra, as well as with the host cluster environment. Our 
BCG sample has the advantage of having detailed X-ray 
values of the host cluster properties available from the lit- 
erature, which BCGs extracted from the SDSS lacks. The 
spectra in this BCG sample also have significantly improved 
signal-to-noise ratios compared to previous samples (Craw- 
ford et al. 1999, and references therein). Most of the previous 
observational studies of BCGs in the UV concentrate on de- 
tecting star formation using the NUV measurements, and 
in mainly cooling flow clusters (Pipino et al. 2009; Donahue 
et al. 2010; Wang et al. 2010; Hicks, Mushotzky & Donahue 
2010, but also see O'Dea et al. 2010 for FUV measurements). 
Here, we investigate the variations in the UV-upturn, which 



in turn can cause scatter in the star formation relations de- 
rived. In Section [2] we describe the UV and optical data. 
In Section [3J we compare the UV-upturn measured for the 
BCGs to that measured for a control sample of ordinary el- 
lipticals matched in mass range and distribution. In Section 
U we test the model predictions by Peng & Nagai (2009), 
which predicts that the UV-upturn phenomenon should be 
most pronounced in BCGs hosted by high mass, dynami- 
cally relaxed, and cool-core clusters. In Section [S] we inves- 
tigate possible correlations between the UV-upturn and the 
luminosity-weighted age and metallicity for the BCG and 
ordinary elliptical samples. We summarise our findings in 
Section [6] 



2 DATA 

We briefly review the relevant selection criteria and proper- 
ties of the resulting sample (detailed in Paper 1) here. The 
sample consists of the dominant galaxies closest to the X- 
ray peaks in the centres of their host clusters. We adopt the 
definition to comply with recent literature (for example De 
Lucia & Blaizot 2007; Von der Linden et al. 2007), where 
the central, dominant galaxy in a cluster is referred to as the 
BCG. According to the above definition, for a small frac- 
tion of clusters the BCG might not strictly be the brightest 
galaxy in the cluster. 

We initially intended to investigate a subsample of 
BCGs with extended haloes (cD galaxies). However, due to 
the difficulties in the classification of cD galaxies and the 
very inhomogeneous definitions in the literature, we cannot 
be confident that all the galaxies in our sample are cD galax- 
ies. Instead, we can say that our sample comprises of the 
dominant galaxies closest to the X-ray peaks in the centres 
of clusters. The sample selection combined three methods 
making the best use of available information from literature 
and astronomical databases: the two well-known galaxy clus- 
ter classification systems that distinguish clusters containing 
a cD galaxy in the centre from other galaxy clusters (Rood 
& Sastry 1971; Bautz & Morgan 1970); an all-sky search 
in the HyperLEDAQ database for galaxies with T-types be- 
tween -3.7 and -4.3 (the de Vaucouleurs Third Reference 
Catalogue classifies cD galaxies as T-type = -4) ; and galax- 
ies that broke the de Vaucouleurs law at large radii. 

In summary: the 49 BCGs in our overall sample were 
classified as cD either in NED and/or have profiles break- 
ing the r 3 law in the external parts. The following global 
criteria were applied to the sample: apparent B-magnitudc 
brighter than 16; distance closer than 340 Mpc; and an ab- 
solute magnitude cutoff at Mb = -20. The sample is not 
complete over any of these parameters, and no UV or X- 
ray criteria were involved in the selection. Nevertheless, the 
galaxies are hosted by clusters with a wide range of X-ray 
luminosities from L x = 0.03 x 10 44 to 17.07 x 10 44 erg s" 1 , 
and it includes BCGs hosted by cooling flow and non-cooling 
flow clusters. 

For this optically-selected sample, high signal-to-noise, 
long-slit spectra were obtained at the Gemini and WHT tele- 
scopes. The spectroscopic data, together with the reduction 

1 http://leda.univ-lyonl.fr/. 
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procedures were presented in Papers 1 and 2, and will not 
be repeated here. The overall BCG sample has 36 galaxies 
with FUV and NUV measurements in the GALEX database 
(ranging from FUV - NUV = 0.20 to 1.32), which will 
be analysed further. This subsample of BCGs (still hosted 
by clusters with a wide range of X-ray luminosities from 
Lx = 0.04 x 10 44 to 17.07 x 10 44 erg s" 1 ) are representa- 
tive of our overall BCG sample. Information from ROSAT 
X-ray data on the properties of the host clusters were col- 
lected from the literature (see SectionQ. Of these 36 BCGs, 
18 have X-ray temperature (Tx) measurements, 27 cluster 
velocity dispersion (<J c i uster ) measurements and 30 X-ray 
offset (Ro//) measurements. All the galaxy and host cluster 
properties used here are listed in Table [ATI 

The GALEX FUV (1344 - 1786 A) and NUV (1771 
- 2831 A) magnitudes and their errors were extracted from 
the GALEX database^ (see Martin et al. 2005). Positional 
cross-correlation of the UV sources in the database was per- 
formed within six arcsec of the optical source coordinates, 
and the UV images were examined by eye to ensure that 
none of the sources were mismatched, and that the sources 
were not near the edge of the detector. In each case, the NUV 
and FUV magnitudes and errors were taken from the longest 
exposures in the GALEX database, and are mainly from the 
All-sky Imaging Survey (AIS). The extracted NUV and FUV 
magnitudes are taken from SExtractor "mag_auto" and are 
the total magnitudes within an elliptical aperture that is 
scaled to 2.5 times the Kron diameter as determined from 
the NUV and FUV images, respectively (Kron 1980; Bertin 
& Arnouts 1996). To minimise aperture effects and instru- 
mental inhomogeneities, we choose to characterise the UV- 
upturn as the FUV-NUV colour, based on GALEX mea- 
surements alone. In addition, there is a lack of enough ho- 
mogeneous V (or r) magnitudes available in the literature 
to keep our sample statistically significant if we characterise 
the UV-upturn as FUV-V. We recalculated the FUV magni- 
tudes so that the effective FUV apertures were equal to the 
NUV apertures, by using the flux measurements from the 
GALEX database and the AB magnitude system of Oke & 
Gunn (1983). The GALEX database defines the NUV and 
FUV magnitudes as mijv = rno — 2.5 log fuv where fuv is 
the count rate and m = 20.082 (NUV) and 18.817 (FUV) 
(Morrissey et al. 2005) . The uncertainty in these zero points 
is ±0.15 mag for both FUV and NUV (Donas et al. 2007). 
We propagate the errors on the magnitudes, extracted from 
GALEX, throughout the study. 

Internal extinction is small in most nearby early-type 
galaxies, and the FUV-NUV colour is essentially reddening 
free for moderate amounts of reddening. We expect a very 
low impact of the internal extinction on the scatter of the 
measurements, and as these corrections are very uncertain, 
we do not include them. 

To quantify the effect of Galactic reddening, the Galac- 
tic extinction was inferred from the 100/xm dust emission 
maps of Schlegel, Finkbeiner & Davis (1998). Using the ex- 
tinction law of Cardelli, Clayton & Mathis (1989), Wyder 

et al. (2007) calculated ^(b-v) = 8 - 24 and e(b-v) = 8 ' 20 > 
which gives corrections in the FUV-NUV colours of <C 0.01 
mag, for the E(B-V) values derived for the galaxies in our 

2 http://galex.stsci.edu/GR4/ 



samplc[f|. Hence, these corrections are negligible and not in- 
cluded here. Similarly, at the low redshift of this sample 
(mean z = 0.0374 ± 0.0029), the if -corrections to the mag- 
nitudes are negligible. 

One of goals of this study is to compare the behaviour 
of the UV-upturn in BCGs with a control sample of ordinary 
elliptical galaxies. Stellar populations in early-type galaxies 
seem to be strongly correlated with the central velocity dis- 
persion, and thus we choose to compare the UV-upturn of 
the BCG sample with a control sample of normal elliptical 
galaxies with the same range of central velocity dispersion 
(which we use as a proxy for mass). For this control sample 
of normal elliptical galaxies, we use the sample of Sanchez- 
Blazquez et al. (2006, hereafter SB06) in the same mass 
range as the sample of BCGs, and excluding any known 
BCGs. The complete SB06 sample consists of 98 galaxies, 
of which 35 belong to the Coma cluster, and the rest are 
galaxies in the field, in groups or in the Virgo cluster. The 
SB06 control sample has 37 galaxies with FUV and NUV 
measurements in the GALEX database. We performed a 
Kolmogorov-Smirnov (K-S) test on the velocity dispersion 
distributions of the two subsamples for which UV measure- 
ments were available within the log a = 2.25 to 2.65 km 
s _1 range (shown in Figure [TJ, where the null hypothesis 
is that the distributions were drawn from an identical par- 
ent population. The two velocity dispersion distributions are 
consistent, and have a test value of 0.260, where a test value 
larger than D — 0.290 indicates that the two samples com- 
pared are significantly different from each other at the 95 
per cent confidence level. We used the same procedure de- 
scribed above to extract FUV-NUV colours for the control 
sample. 

Five galaxies in our BCG sample have GALEX NUV 
and FUV measurements which were used in previous studies 
with similar apertures (Ree et al. 2007; Donas et al. 2007). In 
these studies, the magnitudes were extracted directly from 
the images by the authors. We detect no systematic differ- 
ence between our FUV-NUV measurements and errors and 
those of the previous studies plotted in Figure [3] There is 
one galaxy (UGC00579) in the Donas et al. (2007) sample 
with a very different colour value. We have investigated this 
galaxy without finding any reason for the difference, and we 
have checked that keeping this galaxy included in the sam- 
ple makes no difference to any of the plots and conclusions. 
For the rest of the galaxies the values presented here and 
those measured in the other studies agree within the errors. 

The weak optical emission lines present in some of the 
BCG spectra most likely indicate active star formation re- 
gions or active galactic nuclei (AGN, Edwards et al. 2007; 
Von der Linden et al. 2007) , both of which can alter the UV 
properties of the galaxies. For example, there has been some 
suggestion that relatively tiny amounts of young stars can 
contribute to the observed scatter in the UV-upturn (Han, 
Podsiadlowski & Lynas-Gray 2007; Donas et al. 2007). To 
isolate those galaxies whose UV-spectra may be contam- 
inated by star formation or AGN, we identify the BCGs 
containing emission lines in their optical spectra, using the 



3 When different values such as wff^h =8.24 and = 
8.10 (Atlee et al. 2009) are used, the corrections remain negligible. 
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Figure 1. The velocity dispersion distributions of the SB06 (grey 
- striped) and the BCG (grey - solid) subsamples with UV mea- 
surements. 
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Figure 2. Comparison between the measurements extracted from 
the GALEX archive against previous measurements by different 
authors directly from the images. 



GANDALF routing (Sarzi et al. 2006) as described in Paper 
2. The galaxies with detected emission lines will be indicated 
separately in some of the figures. Of the 36 BCGs studied 
here, eight have been identified as containing emission lines 
in their spectra (see Paper 2). 

We measured the individual emission lines and plotted 
them on a [OIIIJA5007/H/3 against [OIIJA3727/H/3 diagram 
(as [NIIJA6584 and Ha were outside our wavelength range) 
to separate the two major origins of emission: star formation 
and AGN (Baldwin, Phillips & Terlevich 1981; Lamareille 
et al. 2004) in Paper 2. We found all nine galaxies for which 
the H/3, [OII]A3727 and [OIIIJA5007 lines could be measured 
(within at least 2a detections) to be star forming galaxies ac- 
cording to this test. Because the diagnostic diagram used in 
Paper 2 is much less effective at separating different sources 
of ionisation (Stasihska et al. 2006), and because the fraction 
of the current sample containing emission lines is relatively 
small, and with weak emission lines, we do not draw detailed 
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http: / www. leidenuniv.nl / sauron / . 



conclusions about the nature of emission lines in BCGs here. 
It can be seen from Figures [3] [4] [6] and [7] that including or 
excluding the emission line BCGs in this study does not in- 
fluence the conclusions drawn from these plots. It can also 
be seen in Figure [3] that these galaxies are not necessarily 
bluer in FUV-NUV than the quiescent galaxies. 

Several previous studies which observed BCGs in the 
NUV (which is sensitive to the detection of young stars), 
have found a relation between BCGs which show evidence of 
recent star formation episodes and cooling flows in their host 
clusters (Mittaz et al. 2001; Martel et al. 2002; O'Dea et al. 
2004; Hicks & Mushotzky 2005; Rafferty et al. 2008; Pipino 
et al. 2009; Hicks, Mushotzky & Donahue 2010). However, 
not all BCGs in cooling flows show significant star forma- 
tion (Quillen et al. 2008). The UV-upturn is observed in 
the FUV (which is more sensitive to the HB stars), but as 
mentioned above, it is important to assess the influence of 
possible small episodes of recent star formation, as detected 
in the NUV, on our measurements. Rafferty et al. (2008) 
found that active star formation in BCGs depends on the 
presence of cooling flows clusters with short cooling times, 
with the additional requirement of the BCG being located 
close to the X-ray centre (~ 20 kpc). It is therefore also 
of interest to identify the BCGs hosted by cooling flow clus- 
ters, to asses the effects of possible cooling-flow induced star 
formation in the sample. 

We use different symbols in some of the plots to indicate 
galaxies hosted by cooling flow clusters, host clusters with 
definite non-detections of cooling flows, and host clusters for 
which no cooling flow information could be found. Several 
definitions for cooling flow clusters exist in the literature, of- 
ten based on central temperature drop, short cooling time, 
or significant mass deposition rates. We extract the cool- 
ing flow information (see Table [XT]), based on mass deposi- 
tion rates, mainly from White, Jones & Forman (1997), and 
substitute it with information from Rafferty et al. (2006), 
Edwards et al. (2007), and Giovannini, Liuzzo & Giroletti 
(2008). The information on the galaxies in common between 
between these studies agree very well. Of the eight emission 
line galaxies, two are hosted by cooling flow clusters, three 
by non-cooling flow clusters, and three for which the clus- 
ter cooling flow status is not known. Of the 28 quiescent 
BCGs, eight are in cooling flow clusters, 11 in non-cooling 
flow clusters, and nine not known. 

The assessment of the impact of star formation on the 
UV colours is difficult to do, because the main contribu- 
tors to these colours are still a matter of debate. Stellar 
population models including post-asymptotic giant branch 
stars (but not extreme HB stars) indicate that FUV-NUV 
colours are indeed sensitive to the age of the population, but 
not more than FUV-V colours, for example. In a population 
composed of old stars with a small fraction of young stars, 
the FUV-NUV colours are affected only when the age of 
the burst is between 1-2 Gyr (the exact limits depend on 
the fraction of young stars, see figure 12 in Rampazzo et al. 
2007). To conclude, we have tested that even if we take only 
galaxies in non-cooling flow clusters and without emission 
lines, our conclusions hold. 
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Figure 3. The top panel shows the FUV-NUV colour against 
velocity dispersion for BCGs and non-BCGs. The middle panel 
shows the FUV-NUV colour against velocity dispersion for just 
the BCGs, where the filled symbols indicate the quiescent BCGs, 
and the empty symbols indicate the emission-line BCGs. The bot- 
tom panel shows the histogram of the FUV-NUV colour distri- 
butions of the BCGs and non-BCGs. 



3 THE UV-UPTURN AS A FUNCTION OF 
MASS IN BCGS VERSUS NON-BCGS 

Donahue et al. (2010) suggested that their sample of BCGs 
show less dispersion in their NUV-optical colours than that 
of ellipticals as a whole. They posed the question whether 
this uniformity in NUV meant that BCGs are less contami- 
nated from star formation or are more uniform in their age 



and metallicity than their lower mass normal elliptical coun- 
terparts. To answer this question, and since stellar popula- 
tions in early-type galaxies seem to be strongly correlated 
with the central velocity dispersion, as mentioned, we in- 
vestigate the UV-upturn (FUV-NUV colour) of the BCG 
sample relative to the SB06 control sample of ordinary el- 
lipticals in the same mass range. 

We plot the FUV-NUV colours for both the BCG and 
the control samples in the top panel of Figure [3] and the 
histograms of the distributions in the bottom panel of Fig- 
ure [3] Straight line fits to the velocity dispersion plots were 
made for both samples with a least-squares fitting routine, 
and are also shown in the top panel of Figure [3] Statisti- 
cal t-tests were run to explore the presence of correlations. 
Even though neither of the two samples form a statisti- 
cally significant correlation with velocity dispersion (mass) 
at a 95 per cent confidence level, the BCG sample are much 
less scattered with velocity dispersion than the non-BCGs 
over the same mass range. The standard deviation on the 
mean colour values of the BCGs is 0.33 compared to 0.49 
for the SB06 samplH The middle panel shows the FUV- 
NUV colour against velocity dispersion for just the BCGs, 
where the filled symbols indicate the quiescent BCGs, and 
the empty symbols indicate the emission-line BCGs. The 
outliers in this plot are discussed in Section [S] We do not find 
a clear distinction between the active and quiescent BCGs, 
which is in agreement with the findings of Burstein et al. 
(1988) for normal elliptical galaxies. 

To eliminate the possibility that the absence of cor- 
relations is due to the small range in velocity dispersions 
(masses), we used the FUV-NUV measurements of ellipti- 
cals presented by Donas et al. (2007). They found a clear 
correlation over a larger velocity dispersion range and used 
the same apertures as used here. We selected all the ellipti- 
cals in their sample (which includes some ellipticals in the 
centres of clusters) in the same velocity dispersion range 
used here (log a — 2.30 to 2.65 km s _1 ) and fitted a correla- 
tion and t-test as above. In this limited velocity dispersion 
range, we still find a correlation between FUV-NUV colour 
and velocity dispersion for their sample. 

The histogram in Figure [3] indicates another difference 
in the distributions of the FUV-NUV colours of the BCG 
and non-BCG samples, with the ordinary ellipticals having 
larger FUV-NUV colours (i.e. having a smaller UV-upturn). 
We performed a K-S test on the FUV-NUV colour mea- 
surement distributions of the two samples, where the null 
hypothesis is that the distributions were drawn from an iden- 
tical parent population. We find that the distributions of the 
two samples compared are significantly different from each 
other at the 95 per cent confidence level (this has also been 
confirmed with a chi-squared test). The mean FUV-NUV 
for BCGs is 0.791±0.055, whereas it is 1.062±0.080 for the 
ellipticals. 

To further test the validity of these findings, we also 
compare the UV colours of the BCGs with the UV colours 
of the Donas et al. (2007) ellipticals in the same apertures. 
We also plot the 49 ellipticals with logo - > 2.30 km s _1 
from the Donas et al. (2007) sample in the top and bottom 



5 This is the intrinsic scatter since the standard deviations ex- 
pected from the mean errors on the colour values are negligible. 
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panels Figure [3] The UV colours of the Donas et al. (2007) 
ellipticals in the same mass range are also more scattered 
compared to the BCGs (with a standard deviation on the 
mean colour values of 0.42 compared to 0.33 for the BCGs). 
We again find that the ellipticals have bigger FUV-NUV 
colours than the BCGs (we find a mean FUV-NUV colour 
of 1.089±0.059, compared to 0.791±0.055 for the BCGs, and 
to 1.062±0.080 for the SB06 ellipticals). We find that the 
UV colour distributions of the two elliptical samples agree 
very well, whereas the distribution of the BCG sample is 
significantly different from both elliptical samples. 

It is known that BCGs have larger sizes and lower op- 
tical surface brightness than non-BCGs of the same mass 
(Von der Linden et al. 2007). Since we are measuring the 
UV fluxes in a radius scaled to the Kron radius (which itself 
is roughly equal to the half-light radius) , it could be that we 
are measuring a slightly bigger FUV and NUV aperture for 
the BCGs than for the non-BCGs, which could influence the 
NUV-FUV colour measurements. Wang et al. (20f0) simu- 
lated the differences arising in the NUV-r measurements of 
BCGs purely because of their structural differences. They 
found the difference in the NUV-r measurements (between 
the inner, R50, and outer, Rgo, galaxy apertures) to be O.f 
mag bigger for the BCGs than the non-BCGs, purely as a 
result of their structure. The corresponding difference in the 
FUV-NUV measurements can contribute to our result that 
the BCGs have a higher UV-upturn than the non-BCGs, 
but it can not fully explain the large differences in the dis- 
tributions. 

Wang et al. (2010) reported no difference between the 
FUV-NUV colours of BCGs and non-BCGs from SDSS 
data, however they cautioned that this might be ascribed to 
the lower quality of the FUV images used. It has to be noted 
that optically selected BCG samples have shown slightly dif- 
ferent properties from X-ray selected samples (see Burke, 
Collins & Mann 2000; Wang et al. 2010), and the properties 
derived from the two different samples can usually be recon- 
ciled if luminosity-dependent evolution is taken into account. 
However, this does not explain why our result differs from 
Wang et al. (2010), as their FUV-NUV colour comparison 
was also based on an optically selected sample, and both 
samples define BCGs as the brightest galaxy closest to the 
centre of the cluster potential. Wang et al. (2010) also finds 
this null result for both their inner, R50, and outer, R90, 
galaxy apertures, thus the difference can not be as a result 
of different apertures. 

We have also investigated how the scatter and range 
of the FUV-NUV colours in our sample compare to sam- 
ples of BCGs with definite recent star formation. Most of 
the previous samples, such as the cooling flow hosted BCGs 
which showed evidence of recent star formation in Pipino et 
al. (2009), are at intermediate redshifts (z ~ 0.2 to 0.5) and 
should not be directly compared. We identify the nearby 
BCGs (z < 0.1) with blue cores in the optical data pre- 
sented by Rafferty et al. (2008) , and find five of these BCGs 
with reliable FUV-NUV measurements (calculated for the 
same apertures as used here). The FUV-NUV colours of 
these five galaxies range from -1.01 to 0.51 mag (average 
0.06, with a standard deviation of 0.62). If we compare this 
value with the average (and standard deviation) found for 
the FUV-NUV colours of cooling flow hosted BCGs in our 
sample (0.68 ± 0.28), and for non-cooling flow hosted BCGs 



(0.80 ± 0.31), then we do not seem to see a strong influence 
of possible star formation connected to cooling flow clusters 
in our FUV-NUV measurements. Hence, the relative homo- 
geneity of the FUV-NUV colours in the sample here might 
be used as evidence against sporadic episodes of star for- 
mation, since the range and scatter of FUV-NUV colours 
expected for such episodes is large compared to the scatter 
in our BCG sample. We do not analyse this further as our 
main focus is the variations in the UV-upturn (caused by 
HB stars). 



4 ORIGIN OF THE UV-UPTURN: TESTING 
MODEL PREDICTIONS 

Stars with high helium abundance can become UV-bright 
more easily, and therefore, it is believed that helium abun- 
dance can make the UV-upturn more pronounced (Dorman, 
O'Connell & Rood 1995; Yi, Demarque & Oemler 1997). 

It has been suggested that helium sedimentation oc- 
curs in the intracluster medium (ICM) of galaxy clusters 
(Abramopoulos, Chanan & Ku 1981; Gilfanov & Syunyaev 
1984; Qin & Wu 2000; Chuzhoy & Nusser 2003; Ettori & 
Fabian 2006; Peng & Nagai 2009). Under the influence of 
gravity, the heavier helium nuclei accumulate at the cen- 
tres of massive galaxy clusters and can produce significant 
amounts of He abundance (AY > 0.25) in the central re- 
gions of hot, massive clusters. If the UV-upturn is produced 
by extreme HB stars, the helium sedimentation process may 
help to make it more pronounced in those galaxies in the 
centres of the clusters. The UV-upturn phenomena should 
be most pronounced in high mass (or X-ray temperature), 
low redshift, and dynamically relaxed systems according to 
the Peng & Nagai (2009) model. The effect of helium sed- 
imentation, and thus UV-upturn, should also be larger for 
BCGs in cool-core clusters than in non-cool core clusters. 

In particular, the Peng & Nagai (2009) model predicts: 

(i) The UV-upturn correlates with X-ray temperature 
(Tx), and the velocity dispersion of the host clusters 

(o cluster ) ■ 

We use Tx and o c i uster as proxies for the mass of the 
host cluster, and use the values from Table |A"T1 to plot Fig- 
ure U All the values are from spectra observed in the 0.1 
- 2.4 keV band with ROSAT, and using the same cosmol- 
ogy, namely the Einstein de-Sitter model of Ho = 50 km 
s _1 Mpc -1 , Q m = 1 and Qa = 0. We do not convert to the 
Concordance model, as the X-ray temperature and cluster 
velocity dispersion is independent of the cosmological model 
assumed (White et al. 1997; Bohringer et al. 2004). All the 
values for X-ray temperature (Tx in keV) are from White 
et al. (1997). We find no clear evidence that the FUV-NUV 
colours correlate with Tx or o duster in Figure [3] 

(ii) The UV-upturn is more pronounced in BCGs than in 
non-BCGs, and in BCGs hosted by cooling flow clusters than 
hosted by non-cooling flow clusters. 

The former was plotted in Figure [31 and the latter is plot- 
ted in Figure[5] Figure [3] showed significant differences in the 
distributions of the FUV-NUV colours of the BCG and non- 
BCG samples, with the non-BCGs having more pronounced 
differences in FUV-NUV (i.e. a lower UV-upturn and much 
more scatter, as discussed in Section [3}. We find no differ- 
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ence between the FUV-NUV colours of the BCGs in cooling 
versus non-cooling flow clusters in Figure [5] 

It is interesting to note that Wang et al. (20f 0) found that 
among 21 X-ray selected BCGs (12 in cooling flows), those 
located in clusters with central cooling times of less than 
1 Gyr are significantly bluer than those located in clusters 
where the central gas cooling times are long. However, this 
was based on NUV measurements. Unfortunately, in most 
cases we lack the necessary X-ray information to determine 
the cooling time. 

(iii) The UV-upturn correlates with the dynamical state 
of the clusters. 

This is predicted as major mergers, for example, can de- 
stroy helium-rich cluster cores, producing some scatter in 
the UV-upturn among BCGs. Numerical simulations pre- 
dict that the offset of the BCG from the peak of the cluster 
X-ray emission is an indication of how close the cluster is to 
the dynamical equilibrium state, and that this decreases as 
the cluster evolves (Katayama et al. 2003). Thus, the FUV- 
NUV colours were plotted against the offset between the 
X-ray peak of the cluster and the BCG (R g) in Figure [5J 
where R h is in Mpc and from Table [SI] We find no evi- 
dence of a correlation between X-ray offset and FUV-NUV 
colour. 

The host clusters of our BCG sample cover a wide range 
in the parameters chosen to characterise the host clusters. 
The masses of the host clusters range from log a c i ust er = 2.3 
to 3.1 km s _1 , which includes groups to massive clusters. 
Similarly, the range of Tx and 7? Q ff is sufficient to detect 
possible trends (see Rafferty et al. 2008; Bildfell et al. 2008; 
Wang et al. 2010). In all cases, the errors on all the X-ray 
measurements are negligible compared to the scatter. 

We also plotted Figures [4] and [6j distinguishing between 
cooling flow clusters, non-cooling flow clusters and clusters 
for which the information is not known, and we find no cor- 
relations. 

In summary, we do find differences between the FUV- 
NUV colour distributions of normal ellipticals and BCGs in 
the same mass range, as predicted by the models. We do 
not find strong further evidence in our data to support the 
predictions made from the models of Peng & Nagai (2009) 
that the UV-upturn phenomena should be most pronounced 
in BCGs hosted by high mass, dynamically relaxed, cooling 
flow systems as a result of helium sedimentation. We also 
showed in Paper 2 that the BCG sample is not more uni- 
form in age or metallicity than the SB06 sample in the same 
mass range. As such it is still an open question why we 
see differences, such as less scatter, between the FUV-NUV 
colours of the BCGs compared to the normal ellipticals. 

It is also interesting to note that Donahue et al. (2010) 
found no correlation between NUV colour and optical lumi- 
nosity, X-ray temperature, redshift, or offset between X-ray 
centroid and X-ray peak. 

As mentioned in Section [2] we choose to characterise 
the UV-upturn as FUV-NUV colour based on the GALEX 
measurements alone to minimise extinction and aperture ef- 
fects and instrumental inhomogeneities. We note that this 
colour is not optimal for separating the effects of massive 
star formation from those of hot low-mass stars. To test 
whether this choice influences our results, we correlated the 
FUV-NUV colours against FUV-V colours for five BCGs in 
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flow information available in the literature are also indicated. 



our sample which have homogeneous optical measurements 
in Sandage & Visvanathan (1978). We find a strong correla- 
tion between FUV-NUV and FUV-V within the errors (also 
see the correlation in figure 7 of Ree et al. (2007) for their 
nearby galaxies). We find that the scatter introduced in the 
FUV-V values by the large, uncertain aperture corrections 
is larger than the intrinsic scatter in FUV-NUV. 
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5 OTHER CORRELATIONS 

A correlation between the UV-upturn and a property mea- 
surable from integrated light would not only help us to un- 
derstand the origin of the UV-upturn but can also be used 
to utilise the UV-upturn as an age indicator (as discussed 
in Section [T]). 

For normal elliptical galaxies, the UV-upturn is re- 
ported to correlate with Mg, as well as with luminosity 
(Burstein et al. 1988; Donas et al. 2007; Gil de Paz et al. 
2007). In particular, Gil de Paz et al. (2007) found elliptical 
galaxies with brighter Tf-band luminosities (i.e. more mas- 
sive) to be redder in NUV-K but bluer in FUV-NUV than 
less massive ellipticals. 

We plot the FUV-NUV colour against the 2MASS K- 
magnitudes in Figure [71 where we use different symbols to 
indicate emission line BCGs and quiescent BCGs, and BCGs 
hosted by cooling flow and non-cooling clusters, respectively. 
We also plot the 2MASS K-magnitudes of the SB06 sample 
of normal ellipticals in the top panel of Figure [JJ The K- 
band is sensitive to the predominantly red populations in 
massive early-type galaxies and hence is a good measure of 
the total stellar mass. We do not correct the K-magnitudes 
for passive evolution. Using the Bruzual & Chariot (2003) 
stellar population synthesis code with the assumption that 
the galaxies are 10 Gyr old and formed in an instantaneous 
burst, this correction is only -0.2 mag for a galaxy at z ~ 
0.054 (in the jf-band). Thus, it will not make a significant 
difference to whether or not a correlation is found. 

We find no correlation with luminosity, however it is 
interesting to note that when the galaxies with very low 
FUV-NUV colour measurements are excluded (FUV-NUV 
< 0.45 mag), then correlations are found between FUV- 
NUV colour and luminosity as well as FUV-NUV colour 
and mass (Figures [3] and [5]). The low FUV-NUV colour 
measurements are not specific to lower/higher mass galax- 
ies, galaxies with/without emission lines or galaxies hosted 
by cooling/non-cooling flow clusters. It is also not due to ob- 
vious systematic problems in the NUV and FUV measure- 



6 This is clear in the bottom panel of Figure [TJ whereas the top 
panel shows that the BCG outliers are still within the scatter of 
the normal elliptical galaxies. 
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Figure 7. FUV-NUV colour against K-band luminosity. The top 
panel shows the SB06 sample of normal ellipticals for compari- 
son. Clusters with no cooling flow information available in the 
literature are also indicated in the bottom panel. 



ments such as low exposure times for those galaxies. This is 
also the case for the outliers at the high end of FUV-NUV. 

To analyse possible correlations between key absorption 
line indices, such as Mgb, and the FUV-NUV colours, we 
use the indices measured in the central apertures (a c /8jl|. 
We carefully calibrated these measurements to the widely 
used Lick/IDS system of absorption line indices (detailed in 
Paper 2), and plot the H/3 and Mgb indices against FUV- 
NUV colour in Figure [5] We use Mgb instead of Mg2 here, 
since it is less affected by systematic errors for example flux 
calibration, and we re-measure the SB06 indices using the 
same procedure and plot them as reference. As the UV- 
upturn-Mgb correlation is often interpreted as a metallicity 
effect on the UV flux, we also plot the FUV-NUV colours 
against the luminosity- weighted age and metallicity ( [Z /H] ) , 
derived using the predictions of Thomas, Maraston & Ben- 
der (2003) and Thomas, Maraston & Korn (2004). To avoid 
artificial offsets between samples due to the use of different 
techniques to calculate the SSP-parameters, we recalculated 



The 



effective 
(l-e) 



half-light radius was calculated as a c = 
with e the ellipticity (data from NED), r e 



l-c \cos(\PA-MA\)\ 
the radius containing half the light of the galaxy (computed from 
the 2MASS i^-band 20th magnitude arcsec - 2 isophotal radius), 
PA the slit position axis, and MA the major axis. 
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Figure 8. FUV-NUV colour against H/3 and Mgb index measurements, age and metallicity. The green symbols are for the SB06 elliptical 
sample in the same mass range. The blue triangles are BCGs hosted by cooling flow clusters, the red are BCGs hosted by non-cooling 
flow clusters, and the black are BCGs hosted by clusters of which the cooling flow status is not known from the literature. 



these parameters for the sample of elliptical galaxies using 
exactly the same indices and method as used for our sam- 
ple of BCGs. We find no significant correlation between H/3 
or Mgb and FUV-NUV colour (using a f-test), and as in 
Figure [3] we find that the scatter in the FUV-NUV colour 
measurements for BCGs is much less in the plots in Figure[5] 
than for the ordinary ellipticals. As predicted by Yi & Yoon 
(2004) , we do not find a correlation between a-enhancement 
([E/Fe]) and FUV-NUV colour (not shown here). In Pa- 
per 2, we found that BCGs are more metal-rich and have 
higher a— abundance than the control sample of ellipticals, 
and here we have found that the BCGs have a higher UV- 
upturn, however the UV-upturn does not clearly correlate 
with [Z/H] or [E/Fe]. As mentioned in Section [2j exclud- 
ing emission line BCGs from these plots do not affect these 
conclusions. 

Donas et al. (2007) finds a tight anti-correlation be- 
tween FUV-NUV and the Mg2 index for nearby early-type 
galaxies, and Boselli et al. (2005) found a mild trend between 
FUV-NUV and Mg2 in a sample of early-type galaxies in 



the Virgo cluster. On the other hand, Rich et al. (2005) de- 
tected no significant correlation between the FUV-r colour 
and the Mg2 index or the velocity dispersion, using a large 
sample of SDSS early-type galaxies observed with GALEX. 
It is possible that this relation is diluted by lenticular galax- 
ies, as shown in Donas et al. (2007) . This correlation has not 
been previously investigated specifically for BCGs, and even 
though we do not find a statistically significant correlation 
between Mgb (or Mg2 - not shown) and UV-upturn here, 
the properties are less scattered compared to our control 
sample of normal ellipticals. 



6 SUMMARY 

We compare the stellar population properties derived from 
optical spectroscopic data with GALEX UV colour mea- 
surements for 36 nearby BCGs to provide a systematical 
comparison between the UV-upturn in BCGs with that 
of ordinary massive elliptical galaxies. We also test recent 
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model predictions on the origin of the UV-upturn, specifi- 
cally in central cluster galaxies. To do this we analyse pos- 
sible correlations between the FUV-NUV colours and the 
internal properties of the galaxies (velocity dispersion, lumi- 
nosity) and those of the host clusters (mass, distance from 
BCG to X-ray peak, presence of cooling flows). 

We pay particular attention to the model of Peng & Na- 
gai (2009). They pointed out that the helium sedimentation 
process can greatly increase the helium abundance in BCGs 
making them conductive to the formation of UV bright HB 
stars. It has long been suggested that helium sedimentation 
occurs in the ICM of galaxy clusters, and under the influence 
of gravity, the heavier nuclei in the H-He dominated ICM ac- 
cumulate at the centres of the clusters. Thus, measuring the 
UV-upturn in BCGs compared to ordinary massive elliptical 
galaxies should provide a test of this scenario. 

Our main conclusions are the following: 

• We find significant differences between the UV-upturn 
distributions between BCGs and normal ellipticals in the 
same mass range, in that the BCGs have higher UV-upturns, 
and are less scattered than the normal ellipticals. 

• We do not find strong evidence in our data to sup- 
port the predictions made from the models of Peng & Nagai 
(2009) that the UV-upturn phenomena should be most pro- 
nounced in BCGs hosted by high mass, dynamically relaxed, 
cooling flow systems as a result of helium sedimentation. 

• We find no clear correlation between BCG velocity dis- 
persion (mass) and the UV-upturn, but we find the UV- 
upturn of the BCG sample to be much less scattered with 
velocity dispersion than the non-BCGs over the same mass 
range. We also find no clear correlation between UV-upturn 
and luminosity, but the BCGs are less scattered than the 
normal ellipticals. 

• We do not find the UV-upturn-Mg correlation reported 
for some samples of ordinary elliptical galaxies (Burstein et 
al. 1988; Donas et al. 2007), but we do find that the values 
are much less scattered than that of the normal elliptical 
sample that we use as reference. No correlations between 
the age, metallicity and a— abundance of the BCGs with 
UV-upturn were found. 

The UV-upturn is the most variable photometric fea- 
ture associated with old stellar populations, yet it seems 
that the influence of the cluster environment is limited to 
enhancing the UV-upturn in BCGs compared to non-BCGs, 
and independent of the specific properties of the host clus- 
ter (also see Atlee et al. 2009). This absence of correlations 
with host cluster properties makes it difficult to invoke the 
helium sedimentation scenario as the primary mechanism 
responsible for the UV-upturn. 

On both the theoretical and observational grounds, the 
lifetime UV outputs of these old HB stars are very sensi- 
tive to their physical properties (for example helium abun- 
dance). More remarkably, changes of only a few 0.01 solar 
masses in the mean envelope mass of an extreme HB popula- 
tion can significantly affect the UV spectrum of an elliptical 
galaxy (O'Connell 1999). If this interpretation is correct, 
then FUV observations become a uniquely delicate probe 
of the star formation and chemical enrichment histories of 
elliptical galaxies, but only once we understand the basic 
astrophysics of these advanced evolutionary phases so that 
we can explain the diversity in UV observations. 
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Table A2. The velocity dispersions and photometric measure- 
ments of the control sample of ellipticals. 



APPENDIX A: BCG PROPERTIES 

Photometric measurements of the BCGs, and X-ray prop- 
erties and velocity dispersions of the host clusters for all 36 
BCGs are listed in Table I All The velocity dispersions and 
photometric measurements of the SB06 ellipticals used here 
are listed in Table [M\ 
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Table Al. Velocity dispersions and photometric measurements of the BCGs, and X-ray properties and velocity dispersions of the host 
clusters. The BCG a and <r cluster values are in km s _1 and the projected distance between the galaxy and the cluster X-ray peak (R ff ) is 
in Mpc. The * marks at R ff indicate the galaxy is not in the centre of the cluster but closer to a local maximum X-ray density, different 
from the X-ray coordinates given in the literature. The references arc: w = White, Jones & Forman (1997); e = Edwards et al. (2007); g 
= Giovannini, Liuzzo & Giroletti (2008); r = Rafferty et al. (2006); cc = Calculated from Bohringcr ct al. (2000); cl = Calculated from 
Ledlow et al. (2003); cb = Calculated from Bohringer et al. (2004); cr = Calculated from Cruddacc ct al. (2002); st = Struble & Rood 
(1999); s = Struble & Rood (1991); p = Patel ct al. (2006); pe = Perez et al. (1998). All the values for T x arc from White ct al. (1997). 



